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Abstract 
O 

^S| ' Unimodular gravity is based on a modification of the usual Einstein-Hilbert action that allows one to 

(-^ I recover general relativity with a dynamical cosmological constant. It also has the interesting property of 

D ■ providing, as the momentum conjugate to the cosmological constant, an emergent clock variable. In this 

^0 I paper we investigate the cosmological reduction of unimodular gravity, and its quantization within the 

QQ ■ framework of flat homogeneous and isotropic loop quantum cosmology. It is shown that the unimodular 

clock can be used to construct the physical state space, and that the fundamental features of the previous 
models featuring scalar field clocks are reproduced. In particular, the classical singularity is replaced by 
^h| a quantum bounce, which takes place in the same condition as obtained previously. We also find that 

requirement of semiclassicality demands the expectation value of the cosmological constant to be small (in 
{3JQ| Planck units). The relation to spin foam models is also studied, and we show that the use of the unimodular 

time variable leads to a unique vertex expansion. 

(N' 

> 

en 

p 

o 
o 



X 



PACS numbers: 04.60. Pp, 98.80. Cq, 04.60.Ds 



chiou@gravity.psu.edu 



mgeiller@apc.univ-paris7.fr] 



INTRODUCTION 



Loop quantum gravity [l|43| (LQG) is a nonperturbative and background independent quanti- 
zation of general relativity in which the geometry of space-time is treated quantum mechanicahy 
from the very beginning. The philosophy of LQG is based on taking seriously the lessons that we 
have learned from quantum mechanics and general relativity. It relies in particular on background 
independence [j43] and on the connection-dynamics formulation of general relativity [a, loj . The 
application of this theory of quantu m g ravity to symmetry reduced models of general relativity is 
known as loop quantum cosmology |1CH12|] (LQC). In the case of flat A = Friedmann-Lemaitre- 
Robertson- Walker (FLRW) cosmologies with a massless scalar field, the theory strongly departs 
from general relativity in the Planck regime, where quantum gravity effects of LQG become sig- 
nificant and give rise to a new repulsive force [13|. When the matter energy density reaches the 



13,1181, 



20]. The big bang 



and type IX models 



critical density in the deep Planck regime, this force is so strong that the universe undergoes a 
quantum bounce and the big bang singularity is avoided JJL It has been shown recently t hat this 
result is quite robust 15|, and holds in A / cosmologies 16(], A: = 1 closed cosmologies 
k = —1 open cosmologies [l9|, and flat models with an inflationary potential 
sing ularity is also replaced in anisotropic Bianchi type I 2l|, |2j] , type II [2 
24l |. as well as in Gowdy models [25], in which there are infinitely many degrees of freedom. 

In LQC with a massless scalar field (j), the imposition of the quantum Hamiltonian constraint 
on states ^ leads to the evolution equation 9?^ = —0^, where is a second order difference 
operator [unlike in Wheeler-DeWitt (WDW) theory, where it is a differential operator]. There are 
then two ways of constructing the physical state space. The first one is to focus on the constraint 
C = 5^ + = as a whole and use the group averaging procedure 2614281] to map the kinematical 
states onto the physical ones. Alternatively, it is possible to see the scalar field (j) as an internal time 
d la Leibniz in which the evolution of physical quantities are computed in a fully relational manner 
[291 ]. and build a complete set of Dirac observables. However since scalar fields appear in the matter 
Hamiltonian through their momentum squared p^, some difficulties can arise if we consider more 
complicated matter sectors (e.g. with a potential or phenomenological matter), since one needs to 
take the operator square root of the evolution equation in order to see it as a Schrodinger equation 
evolving in internal time (p. In the present work we would like to study the possibility of using 
clocks that are not built with matter fields. We hope that this will help to study the robustness 
of LQC with alternative clock variables, and also shed light on some conceptual issues such as the 
problem of time [30|] . 

All the theories of quantum gravity and quantum cosmology share the same deep philosophical 
questions, the most intriguing and exciting one being certainly the problem of time [31]. By 
this, we refer in a very broad sense to the incompatibilities between the roles played by time in 
quantum mechanics and in general relativity. The notion of time we use in conventional quantum 
theory is deeply grounded in the concept of Newtonian physics, where there is an absolute time 
parameter external to any system, but this conception is simply incompatible with the formulation 
of diffeomorphism-invariant theories. When trying to address the problem of time, it is convenient 
to distinguish among three approaches: those in which a time variable is identified before the 
theory is quantized, those in which a time variable is identified at the quantum level, and finally 
those focusing on the timeless interpretation of quantum gravity. We refer the reader to [3^ for a 
complete review. In this paper we rather focus on an approach known as unimodular gravity, in 
which time emerges at the classical level. 

The idea of unimodular gravity can be traced back to Einstein, who noticed that imposing 
that the so-called unimodular condition det{gf^i^) = — 1 in the action prin ciple leads to a theory 

It has therefore been 
Uii- The 



of gravity in which the cosmological constant remains unspecified ^, l3j] 



argued that this theory might be used to solve the cosmological constant problems 



canonical formulation of unimodular gravity was studied by Henneaux and Teitelboini [39\, and by 
Unruh [40|. They have shown in particular how the cosniological constant arises as a dynamical 
variable accompanied by a canonically conjugate momentum corresponding to a time variable. 



Sorkin has also tried to implement this idea in a path integral approach to quantum gravity 4l|- 
In the generally covariant formulation of Henneaux and Teitelboim, the cosmological constant 
and its conjugate time variable arise through the process of parametrizing [42] a theory with a 
nonvanishing Hamiltonian. It is this point of view that we are going to adopt in order to investigate 
how the emergent time variable can be used as an alternative time variable in LQC. General issues 
concerning the use of the unimodular clock at the quantum level have been discussed by Kuchaf 



43] , and Unruh and Wald [4j] . 



The aim of this paper is to investigate whether the predictions of LQC can be reproduced when 
a clock variable other than a scalar field is used to construct the physical sector and compute the 
evolution of observables. Additionally, we would like to see how the choice of the clock variable 
influences the sum over histories framework that reproduces in spirit the vertex expansion of spin 
foam models. In section |TI1 we briefly review the covariant formulation of unimodular gravity and 
introduce the unimodular time variable. In section IIIIl we investigate its cosmological reduction 
in the connection-triad variables and construct the classical theory with different matter sectors. 
Section [IV] and |V] are devoted to reviewing the quantization procedure of LQC, and the usual 
construction of the physical Hilbert space with a scalar field clock. In section IVH we use the 
unimodular clock variable to study the quantum dynamics of simplified LQC in the pure (without 
matter) de Sitter case. The issues of singularity resolution and occurrence of the quantum bounce 
are addressed in section IVIH where the theory with a nontrivial matter sector is studied. The 
relation with the vertex expansion of spin foam models is discussed in section IVIIII We conclude 
and suggest some further developments in section IIX[ 



II. THE UNIMODULAR TIME COORDINATE 

We take as our starting point the generally covariant unimodular action of Henneaux and 
Teitelboim given by [39| 



S=— [ d^x 
^1^ Jm 



-g{R-2K + Cm) + 2Kd^T'' 



(2.1) 



Here k = SvrG, Cm is the Lagrangian for the matter sources, A a scalar field that serves as a La- 
grange multiplier, and r'^ = {t^,t'^) a space-time vector density. Now consider a fixed coordinate 
system called unimodular coordinates, in which the volume four-form is independent of the four 
metric g'^j^, along with a second coordinate system in which the volume four- form d'^x \/—g appear- 



ing in the action (12. ip is calculated. It has been shown by Kuchaf 43f] that the divergence dyj^ of 
the space-time vector density is simply the Jacobian that transforms from the general coordinate 
system to the unimodular coordinates. This formulation ensures that the action is invariant under 
the full group of diffeomorphisms. 

Let us look at the equations of motion derived from the action (|2.ip . Performing its variation 
with respect to the multiplier A gives the unimodular condition 

V^ = d^T^. (2.2) 



Varying r^ leads to 9^ A = 0, which indicates that A is a space-time constant that can be identified 
with the cosmological constant. Issues about the sign of A will be discussed in the next section. 
Finally, variations with respect to the metric yield Einstein's field equations with an arbitrary 



cosmological constant. 

The interpretation of the vector density t^ is the fohowing. Let us suppose that the manifold 
M is diffeomorphic to S x R, where S is a compact three-manifold. We can integrate (|2.2p over the 
space-time region IZ d M. bounded by two spacelike surfaces S(ti) and S(i2)j with the parameter 
t G M labeling the foliation of M into Cauchy surfaces. The upshot is 

/ ct^xdi^iT^ = / ct^x \J —g = four volume enclosed between 5](ti) and S(t2)- (2-3) 

Jtz Jn 

As we will see in what follows, the cosmological constant A may be regarded as the momentum 
conjugate to a dynamical variable which may be interpreted as a cosmological time parameter 
associated to a surface S(t). This time coordinate, 

T{t)= [ Sxt\ (2.4) 

is just the four volume preceding the spacelike surface S(i), a function that increases continuously 
along any future directed timelike curve. This shows that the evolution of the vector density r^ is 
a pure gauge, except for the mode zero t^ which satisfies 

r(t2) — T{ti) = four volume enclosed between S(ti) and S(t2)- (2-5) 

The unimodular time coordinate acts in fact as a label of equivalence classes of spacelike hyper- 
surfaces with zero four volume separation, which explains its non local nature. 

III. COSMOLOGICAL REDUCTION 

In this section, we introduce the 3 + 1 description of the unimodular theory of gravity in 
connection-triad variables, and then reduce it to flat homogeneous and isotropic cosmological mod- 
els. We then derive the classical equations of motion with and without matter sectors and discuss 
the presence of singularities. 

A. Hamiltonian framework 

Following the formulation of LQG in terms of the connection-triad variables, we introduce a 
triad e" and its dual cotriad o;^, such that e^Wa = 51- Here Latin letters from the middle of 
the alphabet are Lie algebra indices, Latin letters from the beginning of the alphabet are spatial 
indices, and both are ranging from 1 to 3. The three metric qab on the spacelike hypersurface S is 
related to the cotriad via 

qab = ^'a^i^ij- (3.1) 

The spatial geometry is encoded in the densitized triad 

Ef = \e-''e,,,ujiu^l (3.2) 

which serves as our basic phase space variable along with the sti(2)-valued Ashtekar-Barbero con- 
nection 

< = rj. + 7i^:, (3.3) 



where F^ is the spin connection, X* the extrinsic curvature one-form, and 7 the Barbero-Immirzi 
parameter. With these new variables, the action (j2.ip can be recast in the 3 + 1 form as (451 ] 
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dt I (fx 

s 



—EfAi -Nn- N^Ua - N'Gi - - {NA^ - Af° + r"a„A) 



(3.4) 



where A^ is the lapse function enforcing the scalar constraint 
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(3.5) 



^1] 



where F^^ = daA^ — dhA^ + e^j'^AiA'^ is the curvature of the connection A^, and 2KIK^^ 

K\K^^ — K^Ki. N'^ is the shift vector multiplying the vector (or diffeomorphism) constraint Tia, 
and A^* a multiplier used to enforce the Gauss constraint Qi. We see from this action that r" acts 
as a multiplier enforcing the first class constraint Qa = S^A = 0. The symplectic structure is given 
by 



{Al{x),E^{y)} = K^5l5]6\x,y), {rO(x),A(y)} = K5\x,y). 

The total Hamiltonian is 



H= j Sx 



K 



where the new term in the scalar constraint is given by 

Uk = -A^g. 



(3.6) 



(3.7) 



(3.8) 



It enables us to compute the evolution of our canonical variables. In particular, we get A = and 

^, (3.9) 



." = .^ = A.y5 



^0 



which agrees with our definition of the cosmological time coordinate T{t) given in (j2.4p . 

Now we investigate the symmetry reduced model obtained by restricting ourselves to spatially 
flat, homogeneous and isotropic cosmologies. The spatial manifold S is topologically M^, and since 
it is not compact, we restrict all integrations to a fixed fiducial cell V of fiducial volume Vq. S is 
endowed with a fiducial background fiat metric 

yab = °i^y4^ij, (3.10) 

where °c<;* is a fixed cotriad. Isotropic connections and densitized triads are characterized by 



Ai = V'^'c(t)°^L 



(3.11a) 
(3.11b) 



The three metric is related to the fiducial metric by 

Qab = a' {"qah) = «' (X "<^^o) , (3-12) 

where a{t) is the scale factor, which satisfies \p\ = o^Vq . We also have the relation ^/q = y/°q\p\^''^ . 



Because we are using the fiducial cell V, our phase space coordinates c and p are left unchanged if we 
rescale the fiducial metric °qab- Notice that p can be either positive or negative, which corresponds 
to two orientations of the triad. Since the metric is left unchanged under the transformation 
e" — >• —ef and we do not consider fermions in our model, the physics will also remain unchanged. 
Thus this orientation reversal of the triad represents a gauge transformation that needs to be 
handled. 

Because the model that we describe is spatially flat, the spin connection F^ vanishes and thus 
Al^ = 7-ftr*. Since derivatives of i^* do also vanish, we find that 

2^fa^^ = Y^^'A (3.13) 

and the gravitational part of the Hamiltonian can be expressed as 

Having fixed the gauge and diffeomorphism freedom, the constraints Tia, Qi and Qa are automat- 
ically satisfied. Thus we are left with the sole scalar constraint, and one can show that the 3 + 1 
action (|3.4p rewrites 

S= [ dt(—pc+-Kt + N^J\p\c^--NK\p\^/^^. (3.15) 

J \K7 K K7^ AC J 

The canonical term gives us the Poisson brackets 

{c,rf = ?, {r,A} = Ac, (3.16) 



3 



and the Hamiltonian is 



H= f d?xN( %-^\c^ + -K\p\^lA = f d^xN{n+nA). (3.17) 

We see that the first term is the standard symmetry reduced gravitational part of the Hamiltonian 
constraint, whereas the second one corresponds to matter sourced by a cosmological constant A. 
The interesting fact is precisely that this cosmological constant has the time variable T for conjugate 
momentum, which will enable us to obtain an equation evolving in time T when imposing the 
quantum Hamiltonian constraint on wave functions. 



B. Classical equations of motion 

Given the Hamiltonian (j3.17p . we derive the classical equations of motion for our basic variables. 
It is obvious that A is a constant of motion since its conjugate momentum T does not appear in 
the Hamiltonian. With a choice of lapse A^ = 1, the equation of motion for T is given by 



and that for p is given by 



T = {T,n+'HA} = K^ = \pf/^ (3.18) 



p = {p,n + -Ha} = -^^ = -/He, (3.19) 

3 ac 7 



where the dot denotes derivation with respect to proper time t (the choice N = a = y/\p\ would 
correspond to conformal time). The constraint that the Hamiltonian is vanishing gives us 



<= = ±7\/^, (3.20) 

where the plus and minus signs (it) correspond to the expanding and contracting solutions re- 
spectively as will be seen shortly. Note that, without the inclusion of matter, the Hamiltonian 
constraint demands A to be positivq^. 
Substituting ([3:20]) into ([3J9]l yields 

p = ±2^^\p\, (3.21) 

which gives the solution 

\p\ = \po\exp\±2^rit\ (3.22) 




as a function of the proper time t. As expected, at the classical level, the pure (without matter) 
de Sitter universe described here takes infinite proper time to reach p = or p = oo and does not 
exhibit any singularity. Indeed the matter density remains constant and proportional to A and 
therefore does not blow up at any epoch of the evolution. The same thing happens to the Ricci 
curvature scalar since it is proportional to the cosmological constant. To create a singularity and 
reach it in finite proper time, we have to add to the theory some matter with nonconstant matter 
density. 

On the other hand, combining (|3.18p and (13.2ip yields 

P = ±^=J-, (3.23) 

\/\P\ 

which can be integrated to find the solution of p in terms of T 

|p(r)|3/2 _ |p(ro)|3/2 = ±^/3A(r - To). (3.24) 

The triad p is a monotonic function of T. The solution gives either an expanding or a contracting 
universe according to the plus or minus sign (it) respectively. Note that it takes finite time T to 
reach p = 0, but this does not correspond to a singularity as remarked above. 

Following what is usually done in LQC, we can add a scalar field with Hamiltonian 

2 

'H'P = lj^ + \pf'Ui<f>), (3.25) 

where f7((/>) is a potential, and the Poisson bracket is {</>,p</,} = 1- To describe other matter sources, 
such as radiation, dust or any phenomenological component, it is convenient to consider the general 
constant equation of state Pw = wpw The continuity equation 

pw + ^-{Pw+Pw) = (3.26) 

a 



^ In the anti-de Sitter universe, the space is not globally hyperbolic and we do not have a well-posed initial value 
problem by specifying the initial information on a spacelike slice [46|, |47| . Therefore, the Hamiltonian framework 
with the 3 + 1 splitting is not suitable for the anti-de Sitter solution. On the other hand, with the inclusion of 
matters, both positive and negative values of A are permitted. 
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can be integrated to find the dilution law 

p^ = /3a-3(i+'^) = /?|:p|-3(i+'")/2^ (3.27) 

where /? is a constant factor. Dust has an equation of state lu = 0, radiation is described by 
w = 1/3, a massless scalar field hy w = 1, and finally slow-roll inflation corresponds to w = —1 (as 
well as a cosmological constant term). The Hamiltonian describing this generic case is 

-H^ = a^Pv, = \v?'^Pn, = m'^""'^- (3.28) 

The total matter density is given by 

PK+u, ^PK+Pu, = -+ /3|pr=^(^+'")/^ (3.29) 

K 

which blows up at p = if if > —1 (ordinary matter) and atp = ooiftt;<— 1 (exotic matter). In 
the following, we will show that p = or p = oo corresponds to the singularity for w > —1 and 
u; < — 1, respectively. 



Since the total Hamiltonian % + H.^ + 'H«, is vanishing as a constraint, we have 



Note that in the presence of the matter sector, the Hamiltonian constraint admits both positive 
and negative values of A. The equation of motion for T is again T = |pp", and for the triad we 
find 



p = ±2V^ + ^N-(3«'-i)/2 (3.31a) 



= ±2rwA + ^|p|-{3«,+5)/2, (3_31b) 

For w > —\ (respectively w < —1), when p = (respectively p = oo) is approached, the second 
term inside the square root in (13.31ap dominates the first term, and thus we have 

^ |„|-(3«'-l)/4 
"'~^"V 3 '^' 

which leads to 



p^±2\r—\p\-(^'"-'>'\ (3.32) 

V <J 



|p(i)|3(-+i)/4 _ |p(to)|3(-+i)/4 « ^73:5^ (t - to) . (3.33) 

This shows that it takes finite proper time t to reach p = for w > —1 and p = 00 for w < —1. 
Therefore, with the inclusion of matter, the classical solution encounters the big bang or big crunch 
singularity for w > —1 and the big rip singularity for w < —1. 



On the other hand, in the quantum theory (for both WDW and LQC), choosing w = —2 in 
particular makes the algebra much simpler. In order to compare the quantum evolution with the 



classical one, we consider w = —2, and ()3.31bp yields the solutioio 



I + ^HT^/' - V t + ^\PiTo)\'/' = ±^iT- To) , (3.34) 

where the plus and minus signs (it) correspond to the expanding and contracting solutions respec- 
tively. It follows that 



|p(r)|3/2 - |p(ro)|3/2 = ±^3A + 3/3/.|p(r)|3/2 (r _ To) - ^ (r - ro)^ , (3.35) 

which will be used to compare with the counterparts of the quantum theory. 

IV. QUANTUM KINEMATICS 

In this section we are going to study the loop quantization of the unimodular cosmological 
theory at hand. For this we follow the usual Dirac procedure [48], in which the constraints are first 
quantized and then imposed on the wave functionals. Again, because of homogeneity and isotropy, 
the diffeomorphism and Gauss constraints are satisfied, as well as the constraint Qa- We are left 
with the sole scalar Hamiltonian constraint. 

As in the full theory, there is no operator corresponding the the connection itself [3(]. The only 
well-defined operator is the one corresponding to the holonomy of the connection, and therefore we 
need to go back to the expression (j3.14p of the gravitational Hamiltonian and express the curvature 
and the triad in terms of the holonomies. The holonomy of the connection j4^ along an oriented 
edge A °e° is given by [4^] 

(A) /Ac\ ^ „ . / Xc 



hf = cos ^-j 1 + 2sin ^-j Tu (4.1) 

where 1 is the unit 2x2 matrix and Tk a basis of the Lie algebra su{2) such that 2iTk = cr^, 
where a^ is a Pauli matrix. To impose discreteness in the spirit of the so-called improved dynamics 
scheme of LQC, we choose A to be a specific function p,{p) of the triad given by [13l | 

jl^\p\ = A = 2V3TT-flli, (4.2) 

where /pj = Gh is the Planck length squared, and A^ = A the area gap of LQG corresponding to 
the smallest eigenvalue of the surface operator. For the sake of convenience, we introduce the new 
canonical variables 

&^-^, ^ = sgn(p)i^, (4.3) 

which have canonical Poisson bracket 

{&,^} = ^, (4.4) 

and satisfy fic = Xb. 

Now we choose to work in the v representation, where states are given by wave functions ^{i^)- 
The kinematical Hilbert space is Sj = i^gj-av <25 i^matt- ^^ ^^^ next three sections, we will study the 



The solution shows that it takes infinite time T to reach p = oo. This nevertheless is a true big rip singularity 
since it is reached in finite proper time t. 



different matter sectors in relation to the clioice of clock variable. For the moment we focus on the 
gravitational sector. The corresponding kinematical Hilbert space is the space of square integrable 
functions on the Bohr compactification of the real line i^grav = L'^O^Bohi,diJ,Bohr)- Elements of this 
gravitational kinematical Hilbert space are wave functions ^{i') with finite norm 

11^11' = El^('^)l'' (4-5) 

V 

and an orthonormal basis is given by \v) with 

W) = ^vv'. (4.6) 

The operator measuring the volume of the cell V acts like 

Y^{v) = 27r7/|,i|i>|^(i/) = 27r7/f,i|i/|^(i/), (4.7) 

and the holonomies are promoted to operators satisfying 

6^c^(zy) = ^±^^^(i/) = -^{u ± 2A). (4.8) 

Our task is now to promote the Hamiltonian to a quantum operator. To do so, we need to go 
back to the full expression (j3.14p and express the curvature -F^^ and the inverse triad in terms 
of the holonomies. As in the usual improved dynamics quantization of LQC, the action of the 
gravitational Hamiltonian can be expressed as 

7i*(i/, T) = sh^)A{u)sh^)'^{u, T), (4.9) 

where A is a self-adjoint and negative definite operator on S)^^^ defined by 

S7r/2 
A{iy)^{iy) = ^\u\ |i^ + A| - |z^ - A| ^{i^). (4.10) 

In the matter part of the Hamiltonian, the inverse triad |p|~^'^ appears and we need to quantize 
it properly. Thiemann's trick allows us to compute it in terms of holonomies, and its expression in 
the v representation if found to be 
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V-^^{u) = |p|-3/2^(zy) = -^ \v + A|2/3 -\v- A|2/3 ^{v) 

64A3 
= B{u)-^{v). (4.11) 

Let us consider the general case in which a cosmological constant, a scalar field with potential 
[/((?!)), and some matter with equation of state w are present in the matter sector. We just gave a 
prescription to write down the quantum operator corresponding to the gravitational part. For the 
matter part, we are going to use the usual Schrodinger representation, keeping in mind that in the 
unimodular theory, the cosmological constant is also turned into a quantum operator. Requiring 
that the total quantum Hamiltonian operator (for a generic matter sector) annihilates the wave 
functionals ^ E i^^^™ leads to the general equation 

n + ^A^(^)"' + \pIB{^) + U{4>)B{i^)-' + m^T^ ^ = 0. (4.12) 

Now for the sake of simplicity, we shall adopt the two simplifications of the Hamiltonian constraint 

10 



which render the model exactly soluble and simplifies the derivation of the WDW limit 15| . This 
can be done by using the following simpler expressions for the operators A(i') and B{u): 

iM._5^M, i^M = ^. (4.13) 

The first simplification is an extremely good approximation, and it is exact if the wave function 
^{v) has support in the superselected sector v = 4n\/A^pi, n G Z. The second approximation is 
a bit more restrictive, but it has been shown that the error decreases extremely rapidly for high 
values of v/\ [l5( . 

Finally, we require that physical states lie in one of the irreducible representations (either 
symmetric or antisymmetric) of the gauge transformation corresponding to the orientation reversal 
of the triad 11 : v i — t- Ti{v) = —v. Since we do not consider fermions in the present model, wave 
functions are assumed to be symmetric: Yi^{v) = ^{—i/) = ^{i'). 



V. SCALAR FIELD CLOCK 

In this part we investigate the construction of the physical sector of LQC. To begin with, we 
forget about unimodular gravity and review the procedure that is usually followed when matter 
clocks such as scalar fields are used to describe evolution |13|. 

In standard treatments of LQC, the matter sector corresponds to a free scalar field, which is 
monotonic on all classical solutions (this is also the case when k = 1 and A 7^ 0), as can be seen 
from the classical Hamiltonian equations of motion. This suggests that the scalar field is well 
suited to be an internal clock with respect to which we can evolve the other physical variables. 
If the scalar field has a potential, we have p^ 7^ 0, and its evolution does not remain monotonic. 
Thus we have to reset the scalar field clock at different instants of parameter time if we still want 
to use it as a clock. The free scalar field is quantized in the usual Schrodinger representation on 
Sj^^'^ = L'^(R,d(l)) given by 

0^(0)=</.^(</.), p^xI/(,^) = -in^|M. (5.1) 



If we choose a gauge in which the lapse A^ = 1 and adopt the ordering used in [15| , the action of 



the quantum Hamiltonian operator on wave functionals ^(j^, </>) gives the evolution equation 

-^ = -3-G.^.^^^(.,^) 

= ^i^[(i^ + 2X)^{u + 4A) - 2iy^{i^) + (z^ ^ 2A)*(i/ - 4A)] 
4A 

= -Q^{iy,^). (5.2) 

G is a positive self-adjoint operator on i^g^av = -^^(I^Bohr, B{i')dfiBohv)- For each e € [0, 4A), there is 
a Hilbert space ^^™ C i^*^™ in which states have support on the lattice points i^ = e + 4nA, n G Z. 
Since each of these subspaces is preserved under the evolution, there is super selection, and in the 
remainder of this work we will focus on the e = lattice. The difference operator Q acts only on 
the v argument of ^{iy,(j)), and the form of the constraint equation indicates that it is natural to 
view 4> as an emergent time. Elements ^{v, (p) G S)^^^ = i^grav '^■^^™ ^^^ wave functions with finite 

11 



kinematical norm 

ll^lP= f d<P^\^ii^,^)\^ (5.3) 

and an orthonormal basis in S) is given by | u, (p) with 

{iy,^\u',(l)')=6,,,6{(P,(l)'). (5.4) 

Now we can build a set of complete Dirac observables, the scalar field momentum p^ being a trivial 
one. At a fixed value = 00 of internal time, the volume VI^q is also a Dirac observable. Since the 
space of positive and negative frequency solutions is left invariant under the action of these Dirac 
observables, there is superselection. As in Klein-Gordon theory, we therefore choose to work with 
either set. Then elements ^(z/, 0) of the physical Hilbert space S)^'^^^ are the positive frequency 
solutions to the quantum constraint (|5.2p . i.e. solutions satisfying 



i^^^ = Ven-A). (5.5) 



Notice that this equation is equivalent to (15. 2p because the evolution operator Q is independent of 
(f). Indeed, if we act with the operator d/d(j) on (|5.5p . we obtain 

- i^^^^ = ^ (^^)*(^, 'A) + ^e^(z., 0), (5.6) 

and recover (15. 2p since the first term on the right-hand side does vanish. 
The physical inner product is defined at a fixed instant of time as 

('I'l,^2)phys = - Y^ rT^l(i^,0o)^2(l^,</'o), (5.7) 

and is independent of the choice of (J)q since the operator is self-adjoint and positive definite with 
respect to (. , Ophys given above. Note that the state with support at i/ = yields infinite norm 
and hence is not an element of the physical Hilbert spacqj- Now we explicitly define the action of 
the set of Dirac observables. The scalar field momentum acts like 

p^^iiy, 4>) = -^^^^^ = hVQ^{v, 0), (5.8) 

and the volume at fixed instant of time gives a family of relational observables defined as 

V\^,^{v,<i>) = 2^7/|ie*^('^-'^°)|z>|^(i^,<^o) = 2^7/2je^v^(<^-'^o)|z.|*(i.,,/,o), (5.9) 

which represents the action of the volume operator on the state ^{v,(j)) at (p = (pQ and then its 
evolution from po to (j) with the equation (j5.5p . 



VI. UNIMODULAR CLOCK VARIABLE 

In this section we study the quantization of the pure (without matter) de Sitter cosmological 
model corresponding to the unimodular theory. We take the wave functionals to be of the form 

It should be noted that, as emphasized in |l5|, the decoupling of the state \u = 0) does not imply resolution of the 
singularity. 
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"^{u, T), and the usual Schrodinger representation on -^a ~ L'^{^, dT) is used for the cosmological 
constant part: 

f^>{v, T) = T^>{u, T), A^(z., T) = -nih^^^^. (6.1) 

Using again the simpUfied scheme introduced above, the total constraint {T-i + ?^a)^ = takes the 
form 

. d^{i^,T) 3 si^) sk^CAb) 

- '^^r- = ^n^^x ^^(-'^) (6.2) 

^ [*(j/ + 4A) -2^(i/) + ^(z^-4A)] 



where the difference operator is again self-adjoint and positive definite on i^gj-av- Before studying 
the LQC dynamics of this theory, we would like to first ignore the effects coming from quantum 
geometry and focus on its WDW limit. This will make more transparent the construction of the 
physical state space and the Dirac observables, and facilitate the comparison between the two 
quantizations. 



A. Wheeler-DeWitt limit 

The WDW limit is obtained from the LQC dynamical equation (j6.2p by taking the formal limit 
A — )• 0. In the so-called 6 representation, the constraint equation takes the following Schrodinger- 
like form: 



.9 ^(fe,r) _ 3 ,2. 

dT Kh-f^ 



^^^ = 3-2 ^'^(^ T) ^ m{b, T). (6.3) 



An initial solution at time Tq, satisfies the evolution equation 

^(6, T) = e^^(^-^o)^(6. To). (6.4) 

The space of solutions to (j6.3p can be equipped with the physical inner product 

/oo 
db^i{b,To)Mb,To), (6.5) 

-oo 

which is independent of the fixed instant of time To because the evolution operator is self-adjoint 
with respect to it. Finally, we call the resulting physical Hilbert space Q^ ^^. 

Now let us introduce the set of Dirac observables. The first is the cosmological constant A = 
—inhd/dT. It is a constant of motion. The second is the volume at any fixed instant of time, 
which gives a one-parameter family of relational observables defined as 

V\T,^{b,T) = 27r7/2je^e(T-To)|p|^(5^ j.^)_ (g_g) 

In the b representation, the operator i> acts like —2id/db. For any |^) G ^p^^, we have 

(*,y|To^)phys = (V^|To) 
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d 



4z7r7/|i / d6^(6,ro)^*(6,ro) 



oo 
oo 



d_ 

dh 



e-'^(T-To)^{b,T) 



/oo 
db^{b,To) 
-oo 

(^,y|T^)phys + -(T-ro) / (i6^(6,To)6vI/(6,ro). (6.7) 

T J— oo 



Meanwhile, the expectation value for the cosmological constant is given by 

a^(fe,r) 



/oo 
d6^(6,T) 
-oo 



dT 

oo '^-^ 



/oo 
dh^{b,T)hH{b,T), (6.8) 

-oo 



_3_ 

/ ./ — oo 



which is independent of T as A commutes with 0. Also note that (j6.8p automatically yields 
(A) > (see footnoted]). Furthermore, we also define V A as the square root of A via its spectral 
decomposition, such that V A V A = A. The matrix element of V A is thus given by 

/q />00 

(^,VA^)phy. = (VA) = ^ / db^{b,T)b^{b,T). (6.9) 

( J —oo 



It should be noted that V A ^ v |^|) and (v A ) can be positive or negative depending on the state 
I*). 

Equations (j6.6p and (|6.9p lead to 



{V\t) - {V\to) = V3{Va){T-To). (6.10) 

For a semiclassical quantum state which is sharply peaked (highly coherent), we have the good 



"factorization approximation" (A) ~ (vA) or (vA) ~ ±\/ (A), and thus 



{V\t) - {V\to) « ±V3(^) (T - To). (6.11) 

Therefore, the quantum evolution of the WDW theory closely follows the classical trajectory de- 
scribed by (j3.24p if the state is sharply peaked. However, if the quantum state is not sharply peaked, 
the evolution can be considerably deviated from the classical trajectory (i.e. the backreaction is 
appreciable) . 



B. Loop quantum cosmology 

Now we are ready to investigate the properties of the LQC model with the unimodular clock 
variable. Elements "^{u^T) G S) = iOgJ-^v (SDi^^^" are wave functions with finite kinematical norm 

ll^lP = /"dTV|^(i/,r)p, (6.12) 

Jr 

and an orthonormal basis in Sj^^^ is given by | u, T) with 

{v,T\u',T') = 5,,,5{T,T'). (6.13) 
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Physical states satisfy the evolution equation (j6.2p . and the physical scalar product is given by 

(*l,^2)phys = - V ^l{u,To)^2{l^,To), (6.14) 

vr ^ — ^ 

which is again independent of the fixed instant of time Tg for the reasons mentioned abovq^. A 
complete set of Dirac observables is again given by A and V\to- These two operators preserve the 
solutions to (|6.2p and are self-adjoint with respect to (|6.14p because of the self-adjointness of 0. 



Now let us introduce more rigorously the b representation |l5l. l50l]. In the b representation, the 
state 1^) gives the wave functional 

^ib,T)^m)= ^ (6k)M*)= Yl ^"""^-^{^^T). (6.15) 

v=An\ v=in\ 

as the Fourier transform of ^{v^T). From this we see that ^{b + 7r/A,T) = ^{b,T) is periodic, 
and its inverse Fourier transform is 

^(z.,r) = (z.|^) = - r db{iy\b){b\^) = - r dbe-'''^/^^{b,T). (6.16) 

T^ Jo '^ Jo 

The symmetry requirement \l'(i^, T) = ^(— i^, T) that we have chosen to handle the gauge trans- 
formation n is here equivalent to ^(6, T) = ^(— 6, T). In the b representation, the operator i> acts 
like 

z>*(6,r) = -2i^*(6,T), (6.17) 

the operator sin(A6) acts as a multiplicative operator, and the evolution equation (j6.2p becomes 

.9^(6, T) _ 3 sin2(A6). 



dT Khj^ A2 



-*(5,r) = e^(6,r). (6.18) 



To obtain the physical inner product (I6.14|) in the b representation, we have to use the formula for 
the Fourier transform of a Dirac comb: 

^ E ^f/-^)= E ^-""^'^' (6-19) 

which leads to 



Q,t^^ y Q 



j/=4?iA fc=— 00 



if b,b' £ [0,7r/A). Finally, we obtain 

f-n/X 



^l,^2)phys= / db^i{b,To)^2ib,To). (6.21) 

^0 



■* Unlike the physical inner product defined in (|5.7|) . which excludes the state with support at z^ = 0, according to 
(|6.14l) . neither the state with support at ;/ = nor that with the support at i/ = oo is excluded from the physical 
Hilbert space. This again has nothing to do with the issue of singularity resolution. 
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The matrix elements of the Dirac observables are given by 
(^,y|To^)phys = 4i7r7/2jy" db^{b,To)-^^{b,To) 

r-Tr/A 



(6.22) 



e-i^(T^To)^(^b,T) 



Am^ll.j db^{b,n)-^ 

(*,F|T^)phys + -(T-ro) r db^ib,Tof-^^^cos{Xb)^ib,To), 
7 JO 



A 



and 






A2 



as well as 



(*,A/AM/)phys = — r'\b^>{b,Tf^^^{b,T). 
1 Jo ^ 

Note that, as in the WDW theory, (A) > 0. 

In terms of the expectation values, the evolution is of the type 



(6.23) 



(6.24) 



{v\t)-{v\to) = V3{VaJi 
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A (T - To) 



(6.25) 



Again, if the state is sharply peaked, we can make the following factorization approximation: 

{V\t) - {V\t,) « ±a/3(A) Jl-t^{k) (T - To). (6.26) 



6 



Compared with the classical counterpart (j3.24p , this shows that the expectation value of the cos- 
mological constant has to be small ((A) ^ 7^A^) to avoid strong departure from the classical 
trajectory. In this sense, the requirement of semiclassicality demands a small value of the cosmo- 
logical constant. This is in agreement with the result obtained for the vertex expansion in the spin 
foam formulation of LQC as discussed in 5l| . 



VII. NONTRIVIAL MATTER SECTOR 



In this section we consider the quantization of the unimodular theory with generic matter 
sources. In the simplified LQC scheme and the v representation, the quantum constraint equation 
(I4.12P takes the form 



.d^{v,T) 
' dT 



3 sin(A6) sin^ _ ^^| .|_(.,, _ ^ .,| .|_, _ 1^^\ ^^^^^^^ ^^_^^ 



Kh'-f^ A 
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where 

^'" ^ /i(2^7/2^)-+i' ^ = Sfc^V^- C^-^) 

As our purpose of including matters is to create the classical singularity and see if it is resolved 
at the quantum level, in the following, we will consider the simple case by turning off the scalar 
field and including only one phenomenological matter. 



A. Wheeler-DeWitt limit 

In the formal limit A — )• 0, equation ()7.ip with the scalar field turned off becomes 
M{iy,T) ( 3 . 



QT. y^f^^2^ -cj^|-(-+^))vt(^,r) = e*(^,r). (7.3) 

In the z/ representation, matrix elements of the volume operator are given by 

/oo 
di^^{i^,To)\0\^iu,To) 
-00 

/oo 
dzy'l'(zy,T)e^^(^-^o)|i>|e-*^('^-'^o)*(z.,r). (7.4) 

-oo 

If we use the Baker-Campbell-Hausdorff formula to write 

e^e(T-To)|^|^-.0(T-To) ^ 1^1 + ^(^ _ ^^) [0^ |z>|] - i(r - Tof [G, [9, |i>|]] + . . . , (7.5) 



the development will stop at order n + 1 if the exponent of the operator |i>| in (j7.3p is equal to 
n S N {i.e. if —w G N*). For the sake of convenience, let us choose an exotic form of matter with 
w = —2, and define C(^„=_2) = C = k/S'j/A. Then we find 

3 r^ _ . 3 

(*,y|To*)phys = (*,^|T^)phys--(r-ro) / du ^{u,T)b^{u,T) + -C{T - Tq)'' , (7.6) 

7 J-oc 7 

and for the cosmological constant we get 

(*,AvI/)phys = -i^hj'" dv^{v,T)^^^^ 

Q /"CXD /•CX) 

= — / du^!{v,T)h^^![v,T)-hKC I dv^{v,T)\i)\^!{v,T). (7.7) 

7 J — oo 7— oo 



With the inclusion of matter, (A) can be positive or negative, depending on the quantum state. 
The equations above lead to the evolution 

(V^lr) - {V\n) = V3^V^ + /3kF|t^ (T - Tq) - ^(T - ro)^ (7.8a) 



±^/3./( A + (iKV\T) (T - To) - ^(T - ro)^, (7.8b) 



where the square root of an operator is again defined through its spectral decomposition, and the 
factorization approximation has been used to obtain the second line. The evolution described in 
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(j7.8p closely follows the classical trajectory described by ()3.35p if the quantum state is sharply 
peaked, and thus the big rip singularity is not avoided in the WDW theory. 



B. Loop quantum cosmology 

Without the scalar field, the LQC equation (|7.ip reads as 

_,^!S^= (_±,!=Mi^_cj.r <->)*(., T).e*(..r), (7..) 

and the physical scalar product between two states is given by ()6.14p . For w = —2, matrix elements 
of the Dirac observables are given by 



v=An\ 



jy=4nA 

= (*,f|T*)phys + -(r-To) J2 ^(^,r)^^^y^c^^)^(z.,r) 

--C{T-Tof Y ^(i^,T) [l-2sk^)sh^) ^(i^,T), (7.10) 

and 

(*,A*)phys = -^^f^ E ^(^'^)^%^ (7-11) 

u=An\ 

3 v^ ^, ^, sin(A6) sin(A6) ^ , ^, ^ ^ v^ t/ ^m,,^, ^s 
= — E ^(^^^) \ \ ' ^{iy,T)-hKC Y '^{i^,T)\0\^{u,T). 

u=4n\ u=in\ 

Again, (A) can be positive or negative. 
Prom these two equations we get 



{V\t) - {V\n) = ^/3 ( a/a + PnV\T a/i - ^ (a + PkV\t) ){T- To) 



3Pk f^_f^/^^ ^^y^x I (^ _ ^^^2 (7_i2a) 



PA+w\T 



^^''fl-iP^^±^)(T-Tof, (7.12b) 



4 V Pcrit 

where the total energy density operator at the instant T is defined as 



PA+v.\T = -+mT (7.13) 

K 
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in accordance with the classical counterpart (j3.29p . the critical density is given by 

3 

Pcrit = 2T2' C^-^^) 

and the square root of an operator is defined through its spectral decomposition. For the quantum 
state which is sharply peaked, the quantum evolution closely follows the classical trajectory ()3.24p in 
the classical regime when (pa+«;|t) *C /Ocrit- Towards the (putative) big rip singularity, {V\t) keeps 
growing and so does (/5A+«)|r), taking the evolution away from the classical regime. Eventually, as 
(pA+w|r) approaches Pcr\t-, the factor 



_EM^\^^L_{h±M. (7.15) 

Pcrit / y Pcrit 

in ()7.12bp flips signs and (V^It) gets bounced and starts to decrease. Therefore, the expectation 
value of the total density is prevented from blowing up by the quantum bounce. By the same notion 
of singularity resolution suggested in [19], the classical big rip singularity is said to be resolved and 
replaced by the quantum bounce, which bridges the expanding universe with the contracting one. 
The quantum bounce takes place at the epoch when (/Sa+ioIt) ~ Pcr\ t- T he critical density pcrit 
given in (I7.14P is exactly the same as that for other models of LQC [l3|, ll5|. The fact that the 
occurrence of the quantum bounce is signaled by the total (matter plus cosmological constant) 
energy density also agrees with the LQC model featuring a scalar field [13|. 



VIII. SPIN FOAM FORMULATION 

Spin foam models for LQG were introduced as a sum over histories of quantum spin network 
states in a seminal work by Reisenberger and Rovelli [52l | . Interestingly, many different avenues to 
quantum gravity have converged to spin foam models, which now constitute a very active research 
area 1^, ^3;]. Despite the numerous deep technical results that have been obtained over the years, 
several open issues still remain, and it has recently been argued that LQC can be used to shed 



light on these problems 5l|, |53|. In particular, it has been realized that the Hamiltonian theory 
lying behind LQC can be written as a sum over histories that renders the analogy with spin foam 
models more transparent. Let us first briefly recall the ideas behind this formulation and the key 
results that have been obtained so far. Then we will see how the unimodular time coordinate can 
be used to clarify some issues stemming from the standard LQC vertex expansion. 



A. Vertex expansion with a scalar field clock 

In LQC with a scalar field, there are two distinct but equivalent ways to look at the quantum 
constraint equation (j5.2p and construct the physical state space. The first one is to consider 
the scalar field as an internal clock and then construct the positive frequency solutions to the 
Klein-Gordon equation (|5.5p . In this deparametrized framework, the main object of interest is the 
transition amplitude 

A{vf,(^f-u,,(t),) = {uf\e'^^'^f-'t'^\ui) (8.1) 

for an initial physical state \vi) at internal time (pi to evolve to a final state \vf) at internal time (pf. 
The alternative path is to consider the quantum constraint C = —(9? — as a whole and use the 
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group averaging procedure 26l428l| to map the kinematical states into the kernel of the constraint 
operator. This is more in the spirit of the spin foam approach, where one does not have a preferred 
internal time variable to deparametrize the theory. In this so-called timeless framework, attention 
is focused on kinematical states \vi,(f)i) and \vf,(j)f) in io'^™, and the objects of interest are the 
matrix elements of the group averaged inner product, given by the Green functions 



G{uf,(t)f]Vi,(t)i) 



da{vf,(t)f\ 



AaC 



T^iAi)- 



i.2) 



These two ways of constructing the physical Hilbert space are equivalent, and it has been shown 
in 5l|] that one has the equality G(i//, (/>/; fj, 0j) = A{i'f,(f)f;ui,(pi). The idea of the Ashtekar- 
Campiglia-Henderson expansion is then to use a sum over histories formulation to obtain a series 
expansion for these two quantities which reproduces in spirit the vertex expansion of spin foam 
models. 



The vertex expansions that have been obtained 511] for the transition amplitude and the physical 
inner product are respectively given by 



oo 



A{uf,(l)f;ui,(l)i) = ^ 



M=0 

p 

n 



UM-l,---,!^! 
^m^Vm + 1 



UM^M-l^t^M-l^M-l • • • -"J^li^O ^ 



J-^i {uk - 1)! \dH, 



'-Wf^Wj^ 



nfc-1 P 

E 



gtHjnmWm 



(4>f-<t>i) 



\ rij^mi-^-" 



Hn 



and 



G^[vfA!\n,4>i) = X] 



Af=0 



n 



E e 



^M^M-\^VM-\VM~2 • • • "l^lI'O ^ 



t\{n,-l)l 



dQ 



WkU>k 



nu-1 P 



^ UP 



J-\/&n;mWm(.'t'f-'t>i) 



[ n^^„(e. 



WmWm ^WjWj 



i.3) 



(8.4) 



where H = v0, and the subscript + refers to the positive frequency part of the group averaged 
kinematical states. These expressions of the transition amplitude and the inner product are in 
the form of a discrete sum over paths containing M volume transitions. A key observation is that 
although the objects that were computed are the same, they have two distinct vertex expansions. 
One is indeed featuring matrix elements of the operator y@ whereas the other uses matrix elements 
of 0. This peculiarity can be traced back to the form of the Hamiltonian constraint of LQC and 
the fact that the clock variable that is used for deparametrization is introduced through p?. In 
unimodular theory, the fact that the Hamiltonian constraint is linear in the momentum conjugate 
to the time variable enables us to circumvent this difficulty. To see this, we will focus on the timeless 
framework and reproduce the main steps of the procedure leading to the sum over histories. 



B. Vertex expansion with the unimodular time variable 

In the spin foam approach, one does not have a preferred time variable which serves as a 
relational time, and the emphasis is made on the timeless formalism. We will start here with the 
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full constraint operator 



c = e + i^ = e-A = o, (8.5) 



which is self-adjoint on the kinematical Hilbert space S) , and use the group averaging procedure 
to construct the physical inner product. The procedure is the following. We begin by fixing a dense 
subspace S C Sj , which can be taken to be the Schwartz space of rapidly decreasing functions 
f{v,T). If we denote by efc(z^), k £ (—00,00), orthogonal eigenfunctions of with eigenvalue oo^, 
any /(z/, T) G 5 can be expanded as 

/(^, T) = ^jdkjdK fik, A)e^^^efc(z.). (8.6) 

Distributional solutions to the constraint equation lie in the topological dual S* of S, and the 
inclusion relation between all the spaces involved is 5 C S)^^^ C 5*. The group averaging procedure 



2614281] tells us that the solutions are explicitly given by 

*y(i/,T) = y"(iae'°'^/(i.,T) 

dk I dK5{K-ujk)f{k,K)e'^^ek{v) 

dkj{k,u:k)e'^'^^ek{v), (8.7) 

and they satisfy the evolution equation 

*/(z.,r) = e*®(^-^«)vI/y(i.,To). (8.8) 

Given a distribution ^j G S* associated with f £ S, its linear action on an arbitrary function 
g £ S is given by 



[^f\g) = I dkj dA6{A-ujk)f{k,Arg{k,A) 

dkf{k,ujk)gik,ujk)- (8.9) 



Now given two elements f,g £ S, the scalar product between the corresponding states in 5* can 
be defined as 



{^f,^gUys^i^f\g) = {^g\f). (8.10) 



If we consider the full constraint (|8.5p . the object of interest is this physical inner product. Instead 
of considering general kinematical states f{iJ, T), we are going to focus on the basis of i^'^™. 

Let us choose two states ji^i, Tj) and \uf, Tf) in i^ • The group averaged inner product between 
these two states is given by 

G{vf,Tf-Vi,T,) = fda{ut,Tf\e'''^\i^i,Ti) = j daA{uf,Tf-v,,Ti-a), (8.11) 

where A is the probability amplitude for a kinematical state \vi^Ti) to "evolve" to another state 
\vf,Tj) through the unitary transformation by exp(iQ;C). Since the full constraint C is the sum of 
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two commuting operators acting separately on iDgrav ^^^ -5^™, we can write the amplitude A as 

A{i^f,Tf;i^i,Ti;a) = Agrav(i^/; i^*; a)^T(7/;Tj; a), (8.12) 

where 

^grav(i^/; i^i;a) = {ufle-^'^li^i), (8.13) 

and 

AT{Tf;Ti;a) = {Tf\e'''^\Ti) = f dAe'^'^e'^^^f-^^l (8.14) 

Now we need to write down the expansion of the gravitational amplitude. To do so, we rigorously 
follow the treatment carried out in [5lJ, and express it as a sum over discrete histories. It can be 
shown that the amplitude takes the form 



+ 00 



^grav(^'/;i^i;a) = ^ 



A/=0 



n 



UAf^l,...Ul 
Vm^Vm + l 



I^M'^M-l^'^M-l'^M-2 ■ ■ ■ '^t'lJ'O ^ 



L15) 



d 



L (nfc-1)! ydQ 



E 



WkWk , 



1 lljym \^WmW„i ^WjWj) 



where Ai = 1 is the time interval used for the skeletonization of the Feynman path integral. Now 
if we make the assumption that the sum over M in this expression commutes with the integration 
over a in (jS.lip . we find that the group averaged inner product takes the form 



+ 00 

G{uf,Tf;Ui,Ti) = ^ 
M=0 

P 



E e 

l^m^Vm + l 



UM'^M-l^VM-lVM-2 ■ ■ ■ ^VxUo ^ 



J-i(nt-l 



d 



k 

+ 00 

E 

M=0 
P 



L (nfc-1)! yOQ 



rife — 1 p 



WkWk 



m=l 






E e 

Vm^l^m + l 



UM'^M-l^'^M-l'^M-2 ■ ■ ■ ^l^ll^O ^ 



.ii (n.-l) 



d 



n.fc-l 



fc-i p 



E™ 



g*"i"m™m(2/ Ti) 



^^-^ [nk Ij- \dQuikWk J rn=l n^ym (©«'m-!«m - ®WjWj) 



il6) 



This expression is to be compared with the expansion (j8.4p that has been obtained with the matter 
clock. 



If one rather focuses on the deparametrized framework and considers states satisfying the evo- 
lution equation 



^?^^Bn.n 



il7) 
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it can be shown that the vertex expansion of the transition amphtude 

A{uf,Tf;ui,Ti) = (i//|e^®(^/-^'Vi) (8.18) 

is exactly the same as (j8.16p at any order. Therefore not only do we have A{vf,Tf\Vi,Ti) = 
G{vf,Tf;i>i,Ti), but now the two vertex expansions do coincide. This can be traced back to the 
fact that in the unimodular theory the time variable is introduced via a linear function of its 
conjugate momentum. Thus, there is one unique vertex expansion and its expression is simplified 
since there is no square root involved. 

IX. CONCLUSION AND DISCUSSION 

At the classical level, unimodular gravity is equivalent to Einstein's general relativity. The only 
difference resides in the fact that it allows for a dynamical cosmological constant. The price to pay 
for this particularity is that one has to extend the phase space in order to add the cosmological 
constant and its canonically conjugate momentum, which turns out to be a geometric time variable 
measuring the four volume between two hyper surf aces. The main aim of this paper was to inves- 
tigate whether it is possible to use the cosmological reduction of unimodular theory as a starting 
point to build the quantum theory in the spirit of LQC. There is hope that this framework might 



be useful to address issues related to the problem of time and the cosmological constant [371], and 
we think that it is important to study simplified models in order to clarify the statements and give 
an idea of how the construction of a quantum theory can be carried out. 

It has been shown that when the cosmological time variable that arises in the unimodular 
theory as space-time four volume is used, the quantum theory appears to be simpler than that 
of usual LQC with a scalar field. By this, we mean that since the Hamiltonian is linear in the 
cosmological constant, which is conjugate to the unimodular time variable, the dynamical equation 
that is obtained by imposing the operatorial constraint at the quantum level is a Schrodinger-like 
equation. To investigate the relevance of this model, we have first studied its classical dynamics 
with a pure cosmological constant, and then with matter with a nontrivial equation of state in 
order to introduce a singularity. In the quantum theory, we have shown that the Wheeler-DeWitt 
dynamics does not lead to singularity resolution. The classical trajectory is followed and the (big 
rip) singularity is not avoided. In the LQC dynamics, however, the singularity is resolved and 
replaced by a quantum bounce that occurs when the expectation value of the total matter density 
reaches the critical value pcnt = 3/(k7^A^). This critical value is in accordance with the one that 



has been previously derived in LQC with a scalar field clock [la . |15[. Moreover, we have shown 
that in order for a classical solution to exist at late times in the framework of LQC, the expectation 
value of the cosmological constant has to be small enough (in Planck units), with an upper bound 
in agreement with what is usually expected from vacuum energy contributions in flat space-time 
quantum field theory. We see this result as a small clue indicating that unimodular theory might 
be useful to study the quantum dynamics of cosmological models. Indeed, it reproduces precisely 
the results that have been obtained previously, and leads to a linear evolution equation, which 
simplifies the study of the dynamics. We have also briefly studied the vertex expansion of LQC in 
the spirit of spin foam models, and shown that it is possible to obtain a unique expansion when 



the clock that is used is not a scalar field. It has been claimed 5l| that different vertex expansions 
can appear whenever there is a degree of freedom in the theory (which has been taken to be a 
scalar field in previous studies) which we can use for deparametrization. In the unimodular theory, 
the timeless and deparametrized frameworks are trivially equivalent because of the linearity of the 
constraint in the momentum conjugate to the time variable. This supports the idea that two vertex 
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expansions will exist if a matter clock is used for deparametrization, regardless of what happens 
in the gravitational sector (e.g. for Bianchi models). 

There are several directions that one could follow from here. First of all, it is possible to 
use the unimodular theory to investigate more complicated cosmological models such as Bianchi 



models, Kasner models or the Schwarzschild interior J54l. l55l|. This might be a way to introduce 
a cosmological constant along with a time variable in order to compute the dynamics and the 
evolution of the variables of interest. Since the scalar field is not used as a clock anymore, it is now 



possible to quantize it in the polymer representation [56|], and thereby obtain a full loop quantum 
description of the dynamics near and through the bounce. Alternatively, if one takes unimodular 
theory seriously, it is possible to further investigate its quantization in the loop representation. 
This approach might help to understand the regularization of the Hamiltonian constraint in the 
canonical approach, and the dynamics in the presence of a cosmological constant. It is indeed quite 
easy to cast the unimodular action of Henneaux and Teitelboim in the Plebanski formulation, which 
can be taken as a starting point for the spin foam quantization. Unfortunately, in four space-time 
dimensions, we do not have for the moment a good understanding of how to build the full quantum 
theory in the presence of a cosmological constant. 
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